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ABSTRACT 

We present a Chandra observation of the candidate BL Lac object EXO 0423.4—0840. The 
X-ray emission from EXO 0423.4—0840 is clearly extended, and is associated with an optical 
early-type galaxy, MCG-01-12-005, at the centre of cluster CIG 0422-09. We do not detect a 
point source which can be associated with a BL Lac, but we found a small radio source in the 
centre of MCG-01-12-005. The cluster gas temperature mapped by the Chandra observation 
drops continuously from 80 kpc towards the centre, and is locally single phase. We measure 
a metallicity profile which declines outwards with a value O.8Z0 in the centre, dropping to 
0.35Zq at larger radius, which we interpret as a superposition of cluster gas and a dense 
interstellar medium (ISM) in the central galaxy. Although the temperature profile suggests 
that conduction is not efficient, the ISM and intra-cluster medium seem not to have mixed. 
The entropy profile declines continuously towards the centre, in agreement with recent results 
on groups and clusters. The radio source appears to have had some effect in terms of gas 
heating, as seen in the small scale (~10 kpc) entropy core, and the asymmetric hard emission 
on the same scale. 

Key words: galaxies: individual: MCG-01-12-005 and EXO 0423.4-0840 — X-ray: galax- 
ies: interstellar medium — radio: galaxies — galaxies: active: BL Lacertae — galax- 
ies:clusters:intergalactic medium:individual: EXO 0422, C1G0422-09 
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1 INTRODUCTION 

EXO 0423.4-0840 was discovered as part of an EXOSAT high 
galactic latitude survey (Giommi et al. 1991). The source was as- 
sociated with the rUpg = 15.9 mag galaxy MCG-01-12-005, a radio 
source (with 6 cm flux density 44 mjy) with no optical emission 
lines. It was classified as a candidate BL Lacertae object on the ba- 
sis of its X-ray luminosity (above 10^^ erg s^^) and its position in 
the aox — otro diagram, and is included as a BL Lac object in the 
Veron-Cetty & Veron (2003) catalogue. Observations with CAN- 
GAROO give an upper limit for the 7-ray flux above 2 TeV of less 
than 4 x 10~^^ photons cm~^ s"^ (Roberts et al. 1998) but the 
results are not conclusive on the nature of the source. Kirhakos & 
Steiner (1990) determined a redshift of 0.0392 for the source they 
designate IRAS 04235-0840B and identify with the HEAO source 
IH 0422-086, which they classify as a type 2 Seyfert, and associate 
with EXO 0423.4-0840. 

MCG-01-12-005 is an early type galaxy close to the cen- 
tre of galaxy cluster CIG 0422-09, a.1 z = 0.039. White, lones 
& Forman (1997) described it as a cooling-flow cluster and give 
a temperature of 2.9 keV (from Edge & Stewart 1991). ASCA 
data analysed by Ikebe et al. (2002) implied a temperature of 
kT = 2.0to'i keV from the Lx ~ T relation, and a 0.1-2.4 keV 
luminosity of (5.3±0.5) x 10*^ erg s^^. The Galactic absorption 
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for EXO 0423.4-0840 is 6.22x10^'^ cm'^ (Dickey & Lockman 
1990), and the TVh reported in X-ray analysis of the galaxy cluster 
is 6.4x10^° cm"^ (e.g Ikebe et al. 2002). 

In this paper we present Chandra observations of 
EXO 0423.4—0840 which give us new insights into the na- 
ture of this source. Throughout this paper, we adopt a redshift for 
EXO 0423.4-0840 of z = 0.039 and a cosmology with Ho = 70 
km s^^ Mpc~^, Qrn ~ 0.3 and Qa =0.7. The angular diameter 
distance to this source is 159 Mpc, and its luminosity distance is 
172 Mpc. At the redshift of the source, 1 arcsec corresponds to 773 
pc. 



2 OBSERVATIONS AND PREPARATION 
2.1 Chandra data 

EXO 0423.4-0840 was observed with Chandra in December 2002 
(Observation ID 03970) as part of a small sample of BL Lac ob- 
jects with relatively weak cores, which included PKS 0521—356 
(Birkinshaw, Worrall, Hardcastle 2002), with the primary objective 
of studying the interaction between the jet and the group or cluster 
X-ray medium. Since the BL Lac core of EXO 0423.4-0840 was 
expected to produce a high X-ray count rate, the observation was 
performed in faint mode and with a small window with the back- 
illuminated chip S3 on the Advanced CCD Imaging Spectrome- 
ter (ACIS). This window was 128 pixels (of 0.492 x 0.492 arcsec^ 
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(a) (b) 

Figure 1. (a): ACIS image of EXO 0423.4—0840 in the energy range 0.4-7.0 keV. The solid and dashed circles indicate the regions used to extract the global 
and background spectra (see text). The rectangular regions were excluded to avoid the CCD boundaiies. The point source detected by WAVDETECT is also 
excluded, (b) X-ray contours from the 0.4-7.0 keV flux image of EXO 0423.4—0840 overlaid on the Digitized Sky Survey (DSS) optical image. X-ray contours 
are drawn at multiples of 1.1 X lO"'^ photons cm~^ s~^ pixel" ^ . The contours were obtained after excising the compact source and the extended structure A 
(Fig.m, and by replacing them with Poisson samples from the local surrounding (see text for details). We note the different X-ray and optical elongations of 
the central region. The contours are flattened to the north-west because of the edge of the ACIS/CCD. The galaxy MCG-01-12-006, is visible to the south-east. 



each) wide, and 1024 pixels long. We processed the level 1 event 
data using CIAO v3.0.2 and calibration database (CALDB) 2.26. 
We removed bad pixels using the supplied bad pixel file, and re- 
tained events of grade 0,2,3,4, and 6. We verified that there are 
no periods of statistically higher background by generating light 
curves in different energy bands, and the useful exposure time is 
14280 s. The online aspect calculator^ shows that no aspect correc- 
tion was necessary for these data, therefore the X-ray, radio and op- 
tical frames should align to the Chandra aspect uncertainty ~ 0.6 
arc sec. 



2.1.1 Spatial analysis 

The analysis of the diffuse emission started with our use of WAVDE- 
TECT to identify possible point sources in the field. The algorithm 
was used with a source detection threshold'^ set to lO"'' (the mini- 
mum value for the algorithm to give an accurate enough detection), 
and only found a sing le point source at RA = 04''25™51''04, Dec 
= -08°33'36'.'89 (J2000). 

The characterisation of the diffuse emission is affected by the 
half-chip observation, which induces border effects to the north- 
east and south-west. We corrected these effects via the exposure 
map, generated using the aspect histogram and the instrument map. 
When appropriate, the counts image was converted into a flux im- 
age by dividing it by the exposure map. 



for the source, which were computed from Chandra calibration 
database 2.26. This release takes into account a correction for the 
degradation of the low-energy quantum efficiency. We extracted 
spectra in successive circular annuli, excluding the point source and 
small scale structures (see Sec. |5}. The background for all annuli 
was chosen to be the same, and was extracted in a circle of radius 
28.5 arcsec (see Fig.[T(a)t centred at RA = 04''25'"42"268, Dec = 
-08°32'35'.'91 (J2000). For the spectral analysis we excluded the 
regions close to the chip boundary. If not otherwise specified, spec- 
tra were fitted with absorbed plasma models in XSPEC vl 1.3. We 
tried fits using the APEC (Smith et al. 2001) and MEKAL (Mewe, 
Gronenschild, van den Oord 1985; Liedahl, Osterheld, Goldstein 
1995) codes; the latter gave better fits. We scaled abundances to the 
standard solar values from Grevesse & Sauval (1998). 



2.2 Radio data 

A-array VLA data for EXO 0423.4-0840 at 1.4, 4.9, and 8.6 GHz 
were extracted from the National Radio Astronomy Observatory 
(NRAO) archive, and processed into maps through the normal pro- 
cesses of calibration and CLEANing. At 1.4 and 4.9 GHz it was 
possible to self-calibrate the data to improve the final maps, but the 
high noise and short observation at 8.6 GHz made this impossible. 

All three maps easily detect a bright extended radio source 
near the centre of the X-ray emission. At 1.4 GHz the radio map 
has a synthesized beam of 2.3 x 1.9 arcsec, and does not well de- 
fine the structure of the radio source. However, at this lowest fre- 
quency we also find a diffuse source located near RA = 04''25™55', 
Dec =— 08°34'05" which may (or may not) be associated with 
EXO 0423.4-0840. Unfortunately, this structure lies off the field 
of the Chandra X-ray observation. 

The 4.9-GHz image of EXO 0423.4-0840 (Fig. 2) shows the 
1 Ihttp : //cxc. harvard .edu/ciao/threads/arcseccorrection/^adio emission to be from a small double source. At this resolution 
^ For detlirse ^^^^^^^x^^a^fa^^^d^^^o7dowr^^a^do^de^5it)J>t Ml5iM(^^ thtbpamecBilipoflfiltjtSLflcejwelt igparatedti but it is 



2.1.2 Spectral analysis 

For spectral analysis, we accounted for the position-dependent re- 
distribution matrix file (RMF) and the auxiliary response file (ARE) 
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Figure 2. Image in the energy range 0.4-7.0 keV smoothed by a Gaussian filter of cr = 2 pixels. The insert shows a zoom of the same image with overlaid 
contours of the 4.9 GHz map. The first contour is traced at 0.5 mjy, the last at 10 mjy, and they are logarithmically spaced. The figure shows the complex 
X-ray structure present in the core of the cluster. The radio emission lies in a central deficit in the X-ray surface brightness. P indicates the compact X-ray 
source detected by WAVDETECT. The Arc is marginally significant, while the area labelled with A will be discussed in Sec.lsl 



unclear based on the radio data whether one of the components is a 
core and the other is a one-sided jet or both components are lobes 
symmetrically disposed about the core. 

The 8.4-GHz image has 0.4 x 0.3 arcsec full-width-half- 
maximum (FWHM) but a high noise level, so that we cannot cur- 
rently distinguish which component has the flatter spectrum, or 
whether there is a core present between the two components. Fur- 
ther high-resolution, higher-sensitivity radio data are needed to re- 
solve this issue. 



3 X-RAY MORPHOLOGY 

EXO 0423.4—0840 shows clear diffuse X-ray emission (see 
Fig. [T(a)} . The X-ray peak is at RA = 04''25'"51"0, Dec = 
-08°33'36'.'9 (J2000) which is close (less than 5 arcsec) to the 
NASA/IPAC Extragalactic Database (NED) position of galaxy 
MCG-01- 12-005, and coincident with the point source detected by 
WAVDETECT^. The centre of the cluster of galaxies CIG 0422-09 is 
at J2000 coordinates RA = 04''25'"51"02, Dec = -08°33'38'.'5, 
less than 2 arcsec from the Chandra X-ray peak. The Chan- 
dra image suggests that the bulk of the X-ray emission from 
EXO 0423.4-0840 is not due to the expected BL Lac object, but 
is diffuse in nature and originates from the intergalactic medium 
(IGM) of galaxy MCG-01- 12-005, and intracluster medium (ICM) 
of galaxy cluster CIG 0422-09. 

The optical emission of MCG-01- 12-005 from the Digital Sky 
Survey (DSS; Fig. |l(b)t is elongated in the N-S direction, while the 



X-ray emission in the central regions displays a slight elongation 
in the E-W direction. The X-ray contours become more circular at 
larger distance from the centre. 

A more extended structure lies east of the X-ray peak and 
causes the elongation of the central region of the X-ray contours. 
We applied a Gaussian filter of FWHM 0.98 arcsec to the (0.4-7.0 
keV) image, to enhance the extended structure (see Fig.|2j. Two 
X-ray enhancements to the east of the X-ray peak, labelled with 
Arc, seem to be connected to the compact source (labelled with 
P). However this structure is barely significant. The other extended 
structure (marked as A in the insert of Fig.|2j is significant at 1 Icr 
above the surrounding area. In this Figure we also show the radio 
contours from the 4.9 GHz map. We observe that the radio emission 
is located in a region of lower X-ray surface brightness than its im- 
mediate surroundings. We note that the optical centre of the galaxy 
and the southern peak of the radio emission at 4.9 GHz coincide, 
suggesting that this is a likely location for the X-ray centre of the 
object. Unfortunately, the radio data are not sufficient to define the 
location of the core emission. 

To avoid difficulties arising from the uncertain location of the 
X-ray centre, we excluded structure A and the compact source (P) 
for our analysis of the large scale structure. We will discuss the 
structure at the centre in a separate section. We replaced these two 
structures with their local surrounding mean counts and calculated 
the centroid of the resulting X-ray emission, which was found to 
be at RA = 04''25'"51"26, Dec = -08°33'36'.'9 (J2000). For any 
extended source analysis in the following, this will be fixed as the 
centre of the Chandra data. 



3 At the time of writing NED associates EXO 0423.4-0840 with galaxy 
pair MCG-01-12-006, and gives coordinates for EXO 0423.4-0840 of RA 
= 04''25'"55''56, Dec = -08°34'07'.'5 (J2000). The Chandra data show 
this to be an error, since no X-ray emission is associated with MCG-01-12- 
006. SIMBAD correctly associates EXO 0423.4-0840 with galaxy MCG- 
01-12-005. 



4 CLUSTER EMISSION 
4.1 Surface-brightness profile 

We extracted a surface-brightness profile by integrating photons in 
concentric 1 -arcsec wide annuli (masked by the rectangular regions 
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Figure 3. Surface brightness profile derived by integrating counts in the 0.4- 
2.0 keV energy range. The profile was binned to have a minimum S/N=3 
in each bin. The bottom panels show the residual from the 1-/3 and the 2-/3 
models. Bins 2 and 3 were excluded when fitting. 



in |l(a)t out to 118 arcsec and in the 0.4-2.0 keV energy range. We 
subtracted a local background represented by the (masked) annular 
region between 120 and 130 arcsec. The inner 3 arcsec of the pro- 
file display some fine-scale structure (see Fig.|3j, presumably from 
errors in the interpolation over structures P and A as discussed 
above. The excess emission represented by bins 2 and 3 has little 
effect on the parameters of fits to the overall structure of the cluster, 
but it increases the total ■ We thus excluded these two bins when 
fitting. We binned the surface-brightness profile to have a signal- 
to-noise (S/N) of at least 3a above the background in each bin. We 
then initially fitted the profile with a ID /3-model: 



S{r) = So 



{0- 



-3/3+1/2 



(1) 



where Vc and P are the core radius and the index, and So is 
the central surface brightness. The model in equation Q is not a 
good representation of the data, providing best-fitting parameters 
of So = 2.4 ± 0.1 counts s"^ arcmin"^ (5 = 0.51 ± 0.01, 
rc = 14.0 ± 0.8 arcsec ([10.7±0.6]/ifo^ kpc) for a reduced x^ 
xl = 1.64 = 147.9/90 d.o.f. (the xl rises to 1.95 if bins 2 and 3 are 
taken into account, but [3 and rc do not change within the errors). 



The fit improved when we fitted the profile with a sum of two 
/3 models, 



S{r) 



Si{r) + S2{r). 



(2) 



where Si and S2 are defined as in Eq.Q We required the /3 parame- 
ters of both models to be the same in order to have a smooth power 
law shape for the density profile* at large radii. This model gives an 
excellent fit of xl = 1.03 = 90.2/88 d.o.f., with Soi = 2.25 ± 0.15 
counts s^^ arcmin rd = 10.8 ± 1.8 arcsec ([8.3±1.4] h^g 
kpc), 5*02 = 0.79 ± 0.11 counts s~^ arcmin "^ rc2 = 38.4 ± 6.0 
arcsec ([29.7±4.6]/if(/ kpc), and /3 = 0.71 ± 0.06. For compar- 
ison, the analysis of ROSAT data by Reiprich & Bohringer (2002) 
found a value of 13 = 0.722^;;;^°:^ and rc = Wltfi hjo'^ kpc. 
The /3 values are in good agreement if the two-/3 model is consid- 
ered. In Fig.|3]we show the data points and both the single (3 and 
double (3 best-fit models. The lower two windows show the residu- 
als. We derive the gas density distribution in combination with the 
spectral analysis. 



4.2 Global spectrum 

We extracted the global spectrum of the source by selecting a circu- 
lar region of radius 118 arcsec (see Fig. |l(a)t . To limit spurious ef- 
fects produced by the boundary of the exposed part of the chip, we 
excluded two rectangles of size 260 x 105 arcsec^ centred on RA = 
04''25™53''9, Dec = -08°32'4l'.'8, and RA = 04''25'"48''5, Dec 
= — 08°34!49'.'5 (J2000), respectively with the same orientation as 
the CCD (position angle = 30°). The compact source was also ex- 
cluded. We took background from the usual circle (Sec. 2.1.2; Fig. 

El. 

The source has a 0.4-7.0 keV count rate of 1.27 cts s in 
this region. We binned the spectrum to 60 counts per bin prior to 
background subtraction. The global spectrum (see Fig. |4} shows 
strong emission at the energy corresponding to the Fe L complex. 
Other spectral features are also detected, although some of them 
are instrumental (e.g. the region 1.9-2.1 keV near the Ir edge) and 
are excluded when fitting. We fitted the global spectrum between 
0.4 and 7.0 keV with thermal models in XSPEC (version 11.3.0), 
leaving the fractional abundance relative to solar free to vary, and 
initially adopting a Galactic absorption of A^h = 6.2 x 10^" cm~'^ 
(Dickey & Lockman 1990). 

A single MEKAL model gives a bad fit with x^ = 322.0/148 
d.o.f., a temperature of about 2.8 keV and a metallicity of 1.4 solar, 
with large residuals at low energy. A much improved, though still 
unacceptable fit, with x^ = 234/147 d.o.f., is obtained if the ab- 
sorption is allowed to vary. The resulting A^h= (1.19^fl |^Q) x 10^^ 
cm is roughly double the expected Galactic absorption. 

An alternative explanation for the poor fit at low energy might 
be thought to be systematic errors in the used (version 2.26) calibra- 
tion of the quantum efficiency degradation of the ACIS-S3 chip. To 
test this we repeated the fit but used only data above 1 keV. The re- 
sulting fit (Table 1, column 3) is statistically indistinguishable from 
a fit using the full energy range, and so we concluded that the high 
A^'h is not an artifact of the data processing. 

Fitting a two-MEKAL model with TVh as a free parameter in 
the 0.4-7.0 keV range gives x^ = 211.3/145 d.o.f. The F-test in- 
dicates that the two-temperature model is a better representation of 



* We also tried fits with two separate /3 parameters, but the i^-test indicates 
no significant improvement in the fits when the second /3 is allowed to vary 
independently from the first. 
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Figure 4. Background subtracted global spectrum of EXO 0423.4—0840 
extracted in a circle of 118 arcsec. The folded model is a two-temperature 
MEKAL with absorption as a free parameter. The 1.9-2.1 energy bin was 
ignored when fitting due to uncertainties in the Ir instrumental features. 



the data at 99.9% significance (see Fig.|4j- The best-fitting results 
are listed in TableQ The 0.2-10.0 keV total unabsorbed luminosity 
is 3.7x 10*"^ erg s^^. A model composed of a MEKAL plus a power 
law gives a better fit (x^/d.o.f.=222. 2/145) than the single temper- 
ature representation, but does not improve on the two-temperature 
model (see Table and the spectral index is not well determined. 
We thus do not discuss this model any further. 

The results on the global spectrum suggest the following: 

(i) the gas spectrum is consistent with a two-temperature model. 
Spectral fitting to ASCA data (see e.g. Ikebe et al. 1999, Mak- 
ishima et al. 2001) suggested that a two-temperature model can 
mimic a multi-temperature gas as well as a cooling flow (CF), and 
the contribution of the cD galaxy has to be taken into account. We 
thus fitted the global spectrum with a model composed of a MEKAL 
plus MKCFLOW (Mushotzky & Szymkowiak 1988). The high tem- 
perature of the CF model was fixed to be the same as the MEKAL 
model temperature. The metallicities were forced to be the same 
and the absorbing column density left free. The best-fit has = 
216.3/145 d.o.f. The high temperature is kT = (2.61 ± 0.19) 
keV, the low temperature is poorly constrained, with a best-fit value 
of kT = 0.08 keV (and upper limit of 0.45 keV), the metallicity 
Z/Zq = O.SStoAS' and ^e iVH= (l.nto.H) x 10^^ cm-^ all 
in agreement, at 90% confidence, with the two-temperature model 
described above, though a worse fit. The mass deposition rate is 
10±4 M0 yr-\ 

(ii) the absorption is higher than Galactic by more than a fac- 
tor of 2. As shown later (Sec. 4.3.1) the excess absorption encom- 
passes the whole field of view and not only the cool central part 
of the cluster, suggesting that this excess absorption is of Galac- 
tic rather than intrinsic origin. No CO map in the direction of our 
source is good enough to establish if molecular gas contributes 
the needed column. However, there is a good correlation between 
Galactic HI and the infra-red (IR) emitting dust (see Boulanger et 
al. 1996). Pointecouteau et al. (2004) found a remarkable agree- 
ment between their X-ray fitted A'h profile, and the total A'h pro- 
file derived from this correlation. To test if this interpretation is 
applicable to EXO 0423.4-0840, we looked at the Infra-Red As- 
tronomical Satellite (IRAS) map at 100 fim, in the direction of 
the source (Schlegel, Finkbeiner, & Davis 1998). The map shows 
bright clumpy emission, with strong gradients. The IRAS map is at 
far larger spatial scales than the Chandra field of view, but the 100 
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Figure 5. Top: Projected (grey) and deprojected (black) temperature pro- 
files of EXO 0423.4-0840. The two data sets are shifted for better visi- 
bility. The dash-dotted line is the best fit of Equation|5]to the deprojected 
data points (see Sec. l4.4i : Bottom: Projected (grey) and deprojected (black) 
metallicity profiles. 
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^m flux in the vicinity of the field is ~ 7.2 MJy sr 
relation of Boulanger et al. (1996) then predicts A^'h ~ 1.4 x 10 
HI cm~^. This is in good agreement with our best-fit results from 
the global Chandra spectrum, suggesting the absorption excess is, 
indeed, from molecular gas. 



4.3 Temperature and abundance profiles 

In order to perform spatially resolved spectroscopy, we selected 
concentric annular regions out to 1 18 arcsec, masked by the rectan- 
gles seen in Fig. |l(a)| The first bin is of radius 5 arcsec, just outside 
the bright structures in the centre. Possible variations of the spectral 
modelling within this region are discussed in Section|5| but they do 
not affect our conclusions here. 



4.3.1 Projected quantities 

We obtained a projected temperature profile by fitting each annulus 
with a MEKAL model. We used the same background as for the 
global spectrum. All spectra were binned to 15 counts per bin prior 
to background subtraction, and fitted between 0.4 and 7.0 keV. 
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Table 1. Best-fitting results of the global spectral analysis. 



Parameter 


IT 


IT 


IT 


2T 


IT+IPL 


IT+ICF 




(0.4-7.0 keV) 


(0.4-7.0 keV) 


(1.0-7.0 keV) 


(0.4-7.0 keV) 


(0.4-7.0 keV) 


(0.4-7.0 keV) 


TVh (lO^l cm-2) 


0.62 (fixed) 


1 19+017 


1 20+0-60 


1 42+0.27 
^■^^-0.18 


, ,7+0.49 
'-^'-0.13 


, 77+0.13 
'-■"-0.12 


kTl (keV) 


2.83 


2 3g+o.i4 

^■-'"-0.20 


2 50+0.22 


2 47+0.16 
^■^'-0.17 


2 00+0 - 19 


2 61+0.19 

^-°i-0.19 


Z/Zq 


1.42 


72+0-15 


90+0-24 


91+0.22 
"■^^-0.17 


60+°-3° 

'^-™-0.18 




N* (xlO-3) 


6.12 


„ 45+1.03 


7 77 +0.75 
'■" -0.70 


7 80+0- 70 


7 77 + 1.36 
'■'^-6.80 


7 09+0.69 


kT2 (keV) 














TV* (xlO-3) 








33+0.32 






r (keV) 










1.54« 2.34) 




N (xlO-3) 










0.31(< 1.90) 




kTio^ (keV) 












nns+0-45 

l)-US_o.o8 


N 












10 0+4-0 


x2/d.o.f. 


322.0/148 


234.0/147 


153.5/107 


211.3/145 


222.2/145 


216.3/145 



Abundances are from Grevesse and Sauval (1998). In the two-temperature MEKAL model the abundances of the two 
thermal components were tied. N is the normalisation given by XSPEC, as * K{cm-^) = 10"i*l/(47r(DA x (1 + 
2))2) J risnjjdV, where Dj^ is the angular size distance to the source (cm), is the electron density (cm"^)^ and njj is 
the hydrogen density (cm^^) EiTors are quoted at 90% confidence for one interesting parameter. The energy range used for 
fitting is listed for each fit. 



In order to test if A'^h is constant within the area of the X- 
ray detection, we fitted each annulus with a MEKAL model with 
A^'h, temperature and abundances as free parameters. We found that 
all the column densities are consistent with the value found in the 
global analysis (within la errors). To strengthen this result, we fit- 
ted the spectrum in each annulus by excluding data below 1 keV. 
(This also tests for possible spatial variation of the quantum effi- 
ciency degradation at low energy.) These fits are consistent with 
those in the full energy range. We thus confirm that A'h is a factor 
of two higher than the value predicted by the HI measurements at 
all radii, and that it is constant over a scale ~ 120 arcsec (90 kpc), 
with a global value of Nn = {l.Wto.ll) x 10^^ cm^^ jj^g 
following the column density is fixed to this value. 

We next fitted the spectra again in each annulus, between 0.4 
and 7.0 keV, with an absorbed MEKAL model with A'^h fixed. A 
single-temperature model is an acceptable representation of the 
data in each annulus. However, to test possible contamination from 
a non-thermal source in the inner bin, we added a second compo- 
nent, represented by a power law, when fitting the first bin. The 
additional second component is not required (we also added a sec- 
ond thermal model, and the conclusion is similar). 

The projected radial temperature and metallicity profiles are 
shown in Fig.|5|by the grey data points (crosses). Displayed errors 
are la for one interesting parameter. There is a clear temperature 
drop towards the centre, with a minimum projected kT =1.7 keV. 
The temperature rises almost monotonically with radius to reach a 
maximum of 2.8 keV in the outermost bin, which is in agreement 
with the global temperature found with ROSAT (White et al. 1997). 
This suggests that the Chandra data are giving us a fine-scale mea- 
surement of the cool gas in the core region. 

The projected metallicity profile drops continuously from a 
value of 0.8 Zq in the centre to Z / ZQ=Q.'i5 in the outermost re- 
gion. 

4. 3. 2 Correcting for projection effects 

We accounted for projection effects in the profiles of Fig.|5|by us- 
ing the PROJCT model in XSPEC, and assuming spherical symme- 
try. The volumes of each shell are calculated by the model once the 



inner and outer boundaries of each annulus are added to the FITS 
headers of each spectrum. Since the cluster emission extends well 
beyond our outermost annulus, the emitting volume associated with 
the outermost annulus will be too small. The effect is to underesti- 
mate the brightness of the next annulus in. 

We fitted an absorbed, single-temperature MEKAL model in 
which the temperature and metallicity were left as free parameters 
at each radius. The column density was fixed to the global value as 
before. The best-fitting model gives = 802.6/791 d.o.f. 

The asterisks in Fig. |5| show the deprojected temperature and 
abundance profiles of the cluster. When projection effects are taken 
into account, the temperature in the central bin is significantly 
cooler than its respective projected value, since the hot gas from 
the outer regions seen in projection is accounted for. For the outer 
annuli the projected and deprojected temperatures agree within the 
errors. 

The metallicity declines outwards and within the errors is also 
consistent with a constant value of 0.8 Zq within 40 arcsec and 
then drops to Z/Zq = 0.37 in the outermost bin. 

From the deprojected emission measure it is possible to calcu- 
late the ion and electron density directly from XSPEC. Under the 
assumption that the density is constant within each annulus, the ion 
density is simply: 

n.TiH =N x47v 10" Dl (1 + zf /V (3) 

where TV is the normalisation from the deprojected MEKAL model 
fitting. Da is the angular diameter distance and V is the volume of 
the shell. 

The electron density is then obtained from the relation (see 
e.g. Worrall & Birkinshaw 2004) 



where X is the mass fraction in hydrogen, X = 0.74 for normal 
cosmic (solar) abundances. 

The derivation of the electron density directly from XSPEC is 
a suitable alternative option to the classical method of inverting the 
surface brightness profile into a density profile. The latter method 
requires us to estimate the emissivity function at each radius, and 
this is a strong function of the assumed temperature for a cluster at 
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~ 2.5 keV, even if energies below 2.0 keV are used when generat- 
ing the surface-brightness profile. We will see that the two methods 
are, in fact, in agreement. In the top panel of Fig. |6|we show the 
electron density derived from the best-fit normalisation, assuming 
that the density is constant within each annulus. Since the cluster 
extends well beyond our outermost annulus, we overestimate the 
electron density in the last bin. To correct for this effect we note 
that: 

(i) the density derived from the isothermal /3-model described in 
Sec. 14. H is in agreement with the data points in Fig.|6| for all annuli 
but the last (continuous line in Fig.|6j, despite the clear temperature 
variations (Fig.|5}; 

(ii) the ROSAT data provide the best measures of the cluster at 
larger scale, and the best-fitting j3 found by Reiprich & Bohringer 
(2002) is in excellent agreement with our (i value. This suggests 
that the extrapolation of our best-fitting /3-model is a reasonable 
representation of the outer shells of the gas emission. 

For these reasons, we corrected the density in the last annulus to 
equal the /3-model-derived density at its distance from the centre 
(83 arcsec). We added a 20% systematic error to this bin. The cor- 
rected density is 50% lower than the original value, as shown in 
Fig-IUby the dashed point. A similar argument applies to the pro- 
ton density profile. We used the corrected values to calculate the 
entropy and pressure profiles below. 

We also note that the first data point is higher than (but consis- 
tent with) the /3-model profile (continuous line). This arises from 
the exclusion of the excess emission in the central 5 arcsec when 
fitting the (3 model (second and third bin in Fig.|3}, while this ex- 
cess was not excluded when we extracted the central spectral annu- 
lus (but note that the compact source and the extended structure A 
described in Sec.|3|were masked). 
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4.4 Entropy and pressure profiles 

The middle panel of Fig. |6| shows the profile of 5 = T nZ'^^^ 
(Lloyd-Davies et al. 2000), which is related to the entropy per unit 
mass in the gas, derived from the temperature and density profiles 
obtained from the spectral fits. 

It is useful to fit the temperature profile with an analytical 
function. Here we use (Allen et al. 2001) 



Radius (arcsec) 
1 



T = To + Ti X 



{r/r-c 



(5) 



(1 + (r/rc)") 

where our best-fit has To = 1.22 ± 0.13 keV, Ti = 1.40 ± 0.17 
keV, and rc = (10.86 ± 3.43) arcsec, close to the value of rc from 
the radial profile fit. The exponent, r;=2, was fixed to the best-fitting 
value found by Allen et al. (2001). We then calculated the entropy 
profile using the best-fitting temperature model and the isothermal 
/3-model-derived electron density profile. This is shown in Fig.|6|by 
the dashed line. The dot-dashed line shows the S oc r^'^ behaviour 
found in analytical and semi-analytical models (Tozzi & Norman, 
2001; Voit et al. 2002), and which is also found in numerical simu- 
lations without complex gas physics (e.g Borgani et al. 2001). The 
normalisation is arbitrary. 

A simple power-law fit of the form S = Axr" lo this entropy 
profile yields a slope a = 0.77±0.04 (r is in arcsec). The effect of 
adding a constant (S = A x r" + B) is. to increase the slope. The 
full details of the power law plus constant fit are A = 1.04 ± 0.27, 
a = 0.99 ± 0.07 and B = 3.02 ± 0.73. An f -test indicates that 
there is less than 1 per cent chance that the simple power law is a 
better fit. 
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Figure 6. Top: Electron density derived from the emission measure fitted 
in each annulus, assuming the density to be constant in each bin. The con- 
tinuous line is the density profile derived from the double /3-model fitting of 
the surface brightness profile in the energy range 0.4-2.0 keV. Dot-dashed 
lines are the upper and lower Icr limits of this profile. Middle: Entropy pro- 
file. The continuous line is the best-fitting power law (plus constant) model 
(see text). The dashed line represents the entropy profile derived by using 
the best-fitting temperature profile model and the /3-model-derived density 
profile. The dot-dashed line represents the entropy profile when the expo- 
nent is fixed to be 1.1 as expected from shock heating in spherical collapse 
(e.g. Tozzi & Norman 2001; Voit et al. 2002); Bottom: Pressure profile. The 
continuous line is the pressure profile derived from the deprojected temper- 
ature profile and the density profile implied by the /3-model fit. 
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It is now common to find the entropy profile scaled by r2oo, 
the radius within which the mean cluster density is 200 times the 
mean density of the Universe. If we use the equation in Evrard, 
Metzler & Navarro (1996), and the external temperature of the clus- 
ter (2.9 keV), we obtain r2oo = 1-4 Mpc for EXO 0423.4-0840. 
This implies that our outermost annulus corresponds to ~ 0.05 
r2oo, which is too small to allow useful comparisons with pre- 
vious work on clusters with similar temperatures. Pratt & Ar- 
naud (2004) found 5" oc r"'^''*" " for a small sample of clus- 
ters with kT > 2 keV. Although this is calculated in the range 
(0.05-0.5) r2O0! it looks consistent with the result we find here for 
EXO 0423.4-0840. 

The bottom of Fig. |6| shows the pressure profile, where we 
compare the results when we use the XSPEC proton density profile 
(data points) and the isothermal /3-modeI density profile (continu- 
ous line). 



5 THE CENTRAL REGION 

In previous sections we have described the small-scale structure 
in the central 4 arcsec of EXO 0423.4-0840 (see Fig. |2j. In 
the surface-brightness profile (Sec. l4.lt we found excess emission 
above a /3-model, despite the fact that the compact source P de- 
tected with WAVDETECT, and the more elongated emission labelled 
as A in Fig.|2|were masked before extracting the radial profile. We 
also noted that the centre of the galaxy (and the cluster) displays 
a lower X-ray surface brightness than the immediate surroundings. 
In this section we explore the nature of the compact source and the 
extended structure, and investigate the relation between the central 
gas and the radio source. 

We extracted a spectrum of source P from a circular re- 
gion of radius 1.35 arcsec centred on RA= 04''25'"51''02, Dec = 
-08°33'36'.'90 (J2000). The local background was estimated in the 
surrounding annular region within radii 1.35 and 2.95 arcsec. There 
are only 38 net counts in this region. Adopting C-statistics we found 
that a power law of spectral index 2.5l'lJ | is a reasonable descrip- 
tion of the data (the column density was fixed to the global value). 
With this model the unabsorbed 0.2-10 keV flux of the source is 
3.9x10^^* ergs s^^ cm^^. A thermal model also describes the 
data well, but the temperature is not constrained. With the data 
available, we can not state if the source is a clump of thermal gas 
within the cluster or a background (foreground) quasar, although 
the non-thermal interpretation seems more plausible from the spec- 
tral analysis. Given the positional offset between this source and 
the radio emission, it is unlikely that this is the BL Lac object we 
were looking for. We mask this source for the following analysis. 

We also extracted a spectrum of structure A. We have even 
fewer photons but the spectrum shows a clear thermal signature, 
and a temperature consistent with the temperature that we found for 
nearby gas when this structure was excluded (Sec. l4.3t . Since the 
surface-brightness profile displays excess emission between I and 
4 arcsec from the centre, even after excluding structure A, we sug- 
gest that A is of the same nature as the excess surface brightness, 
and its brightness may be due to projection effects or the particular 
geometry of the source. 

To investigate the spectral signature of the excess sur- 
face brightness, we extracted the surface brightness profile of 
EXO 0423.4—0840 between and 15 arcsec in smaller radial an- 
nuli than in Sec. 14. II and in two energy bands: a soft band (S) at 
0.4-2.0 keV, and a hard band (H) at 2.5-6.5 keV. At larger scale 
the profile was binned in 5 arcsec wide annuli. We then obtained 
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Figure 7. Top: Hardness ratio profile. Radial bins are 0.5 ai'csec wide out 
to 10 arcsec and then 5 arcsec wide. Bottom: Image of EXO 0423.4—0840 
in the energy band 0.4-2.0 keV smoothed with a gauss filter of tr = 0.98 
arcsec. Black contours are obtained from the Gaussian filtered hard image 
with the same cr. White contours are from the 4.9 GHz map. The lower 
contour is traced at 0.5 mJy and they are logaiithmically spaced up to 10 
mjy. 



an hardness ratio (HR) profile by calculating [H — S)/ {H + S) 
(Fig. 0. The average value of the HR is -0.77, reflecting the soft 
spectrum of the gas. Points at radius R < 5 arcsec are consistent 
(within the statistical errors) with this average value. We generated 
S and H images and smoothed them with a Gaussian filter of a=2 
pixels (0.98 arcsec). In the bottom of Fig. we show the soft im- 
age in grey scale. The black contours are obtained from the Gauss 
filtered hard image and the white contours are from the radio map 
at 4.9 GHz. The hard emission is detected above the background 
at 4.5(7. We observe that the hard emission is centred in the region 
of low surface brightness soft X-ray emission coincident with the 
optical centre of the galaxy, and that it is oriented in the same direc- 
tion as the radio source but is more extended. The peaks of the hard 
X-ray and radio emission are not coincident, but this is likely to be 
a consequence of the smoothing, and is not significant. The peak to 
the north-west (black contours) is probably an unrelated source. 

What is the physical relationship between the radio source 
and the surrounding ICM? First of all, there are two peaks in the 
4.9 GHz radio map. We label as S the southern component at RA 
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= 04''25'"51"238, Dec = -08°33.'37'.'51, and as iV the northern 
component at RA = 04''25'"51''325, Dec = -08°33'35'.'96. The 
8.46 GHz map shows a peak of flux density 0.8 mjy between com- 
ponent S and A'^, while components S and are resolved and show 
steeper spectra. More sensitive radio observations would be needed 
to decide whether the core of the source lies in S or A'^, or in the 
weak component that appears to lie between them. If we make the 
assumption that the source is in equipartition between radiating par- 
ticles and magnetic field, we find an equipartition magnetic field of 
7.4 nT and a pressure of 1.5x 10^^^ Pa for both components (we 
assumed a spectral slope of 2.2, minimum electron Lorentz factor 
of 100, maximum electron Lorentz factor of 10'', and filling fac- 
tor equal to 1). From the pressure profile discussed in Sec. l4.4l the 
pressure of the external medium at the distance of each of the ra- 
dio components from the optical centre is ~ 4 x 10^^^ Pa, which 
is only a factor of ~ 2 from the internal minimum pressure of the 
source, suggesting that the source is close to pressure balance. 



6 DISCUSSION 

We summarise below the results from our study of 
EXO 0423.4-0840 using the high resolution of Chandra. 

(i) The X-ray emission from EXO 0423.4-0840 is extended, 
without the strong point-like component expected for a BL Lac 
object. EXO 0423.4—0840 is in fact associated with a cluster of 
galaxies, and an early type galaxy which dominates the central op- 
tical light; 

(ii) we found absorption in excess of Galactic. This is not intrin- 
sic to the source but of Galactic origin, as shown by the agreement 
between our X-ray measure and an independent measure which 
uses the A'^e-IR emission relation; 

(iii) we find that the gas is locally isothermal, with a temperature 
decreasing towards the centre, as found for cooling-flow clusters, 
and with a temperature of ~ 3 keV at r > 40 kpc; however, a CF 
model gives a mass deposition rate of only 10 AIq yr~^; 

(iv) the metallicity profile declines outwards and it is consistent 
with an almost constant Z = 0.8^© at r < SOh^g kpc, and then 
drops to 0.37^0; 

(v) the gas entropy profile behaves as S* oc j-0-99±o o'^^ close 
to the relation predicted by analytical models of shock heating in 
spherical collapse, but shows a small core reflecting possible en- 
tropy modification in the centre; 

(vi) the central radio source has a minimum pressure that is a 
factor of 2 below the external medium. Hard X-ray emission from 
the cluster centre has the same orientation as the radio source, but 
is more spatially extended. 



6.1 ICM or IGM? 

The current Chandra observation gives us a precise look into the 
core of an intermediate temperature galaxy cluster. The radial tem- 
perature profile is reminiscent of cooling core clusters, harbouring 
denser gas in the centre. Are we observing a classical cooling flow? 
Arguments against the classical cooling-flow models have been ex- 
tensively discussed since the advent of XMM-Newton (e.g. Tamura 
et al. 2001; Peterson et al. 2003; Bohringer et al. 2002). In our case, 
fitting a cooling-flow model to the global spectrum does not give a 
better fit than a two-phase model, but the temperature profile does 
suggest a cooling flow. 



Another interesting possibility is that we are observing the in- 
terstellar medium (ISM) of the early type galaxy. Indeed our Chan- 
dra observation maps X-rays only within the optical envelope of the 
galaxy. A clear indication in favour of this interpretation is given by 
the metallicity profile. We are not able to map individual elements 
with the current data, but the fitted metallicity should be dominated 
by the iron abundance. The EXO 0423.4—0840 metallicity profile 
clearly shows a decline from the centre outwards, which can also 
be seen as a break at ~ 35/ifo^ kpc, with an average global rela- 
tive abundance of ~ 0.8 within this radius and Z / Zq=Q.2i1 out- 
side. Work based on ASCA data found that clusters harbouring 
cD galaxies display a significant increase in metal abundance to- 
wards the centre (e.g. Fukuzawa et al. 2000; Matsushita et al. 1998). 
Similar results were found with BeppoSax data (e.g. De Grandi & 
Molendi 2001), and suggest that this is due to metal enrichment 
by SN type la in the cDs (e.g. Finoguenov & Ponman 1999, Mat- 
sushita et al. 2002, 2003). The optical galaxy that we associate with 
EXO 0423.4-0840, MCG-01- 12-005, is not weU studied and it is 
not catalogued as a cD galaxy. However, this interpretation is still 
plausible, since prior to the advent of Chandra it was almost im- 
possible to map in such detail the inner regions of galaxy clusters, 
and objects without the striking optical properties of cD galaxies 
are just starting to be discussed. 

An acceptable estimate of the galaxy X-ray temperature is 
given by the first bin of the temperature profile (kT^ 1.3 keV). We 
find a corresponding X-ray luminosity of ^ 3.6 x 10*'^ erg s~^, 
once projection is taken into account, and after extrapolating to 40 
kpc (when the cluster emission dominates). Although the temper- 
ature is rather high, the values we find are in agreement with the 
Lx — T relation of O' Sullivan et al. (2003) for early type galax- 
ies, adding another piece of support to the ISM interpretation for 
the X-ray emission within 40 kpc of this object. While the temper- 
ature increases rapidly after the first bin, and at 40 kpc it is already 
~ 2.5 keV which is too high for the atmosphere of a galaxy, we 
note that the metallicity and temperature profiles could change in- 
dependently with different radial functions because of the different 
physics associated with the transport of heavy ions outwards and 
thermal energy inwards. 

We also note that the gas is consistent with a two-phase metal- 
licity, with a value close to solar within 35 arcsec and close to a 
common value for galaxy clusters beyond 35 arcsec. This suggests 
that little mixing between the ISM and ICM has occurred. 



6.2 The impact of the radio source 

We can use the comparison between the radio source and the 
ambient medium properties to map the formation history of the 
core of EXO 0423.4-0840. We find that the gas entropy de- 
clines continuously towards the centre, in agreement with recent 
results on galaxy groups and clusters (e.g. Ponman, Sanderson & 
Finoguenov 2003; Pratt & Amaud 2003; Khosroshahi, Jones & 
Ponman 2004, Pratt & Amaud 2004), typically of the same tem- 
perature as EXO 0423.4—0840. Moreover, the large-scale radial 
dependence of the entropy we measure is in agreement with both 
analytical models of spherical collapse, and numerical simulations 
based solely on gravitational collapse (Tozzi & Norman, 2001, Bor- 
gani et al. 2001, Voit et al. 2002). It does seem that there is or 
there has been some modification of the core entropy of this cluster 
(galaxy). This is seen in the radial dependence of the entropy, where 
a simple power-law fit gives a shallower slope than expected from 
pure shock heating, and where a power law plus constant model is 
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clearly a better fit. However there is no flat isentropic core (entropy 
plateau). 

Tiie deviation of other classical scaling laws, such as the 
Lx — T relation (e.g. Arnaud & Evrard 1999), from the predic- 
tion of purely gravitational models has suggested that some extra 
heating has occurred. It has been proposed that sources at the cen- 
tres of clusters (such as AGN or radio galaxies) could cause this 
heating (e.g. Valageas & Silk 1999; Wu, Fabian & Nulsen 2000; 
Churazov et al. 2002; Ruszkowski & Begelman 2002). The small 
core in the entropy distribution of this object would suggest that en- 
tropy has been modified in some way, and the central radio source, 
although weak, may be a plausible candidate for this modification 
(e.g. Roychowdhury et al. 2004). 

An alternative observational indication is given by the com- 
parison between the radio and X-ray luminosity in the cooling re- 
gion, which should be similar, once a reasonable factor for the con- 
version between power and radio emission is assumed. The radio 
source at the centre of EXO 0423.4—0840 is particularly small. 
The radio luminosity is of order 5 x 10''*' erg s~^, while the X- 
ray luminosity of the soft thermal component within 35 arcsec is 
~ 3.6 X 10^^ erg s~^. If only 0.1% of the total power output goes 
into the radio emission, then the rest of the power can heat the gas. 
However, this is a rather low efficiency for conversion of the source 
power into radio emission, and also a high efficiency for heating. 
From this evidence, it is unlikely that the radio source is supplying 
sufficient heat to stop any cooling flow. 



6.3 The age of the radio source 

Since we do not have an independent measure of the magnetic field, 
we assumed equipartition, yielding the result that the minimum in- 
ternal pressure of the source is about a factor of 2 lower that the 
external pressure. The minimum pressure is a lower limit, thus we 
are facing two possibilities for the small size of the source: the radio 
source is old and is dying, or the radio source is young. 

In the first case, the source was bigger and more powerful in 
the past. It may have heated the external medium which is now 
hotter than at the origin, as can be inferred from the relative distri- 
bution of hard and soft emission in the core of EXO 0423.4—0840 
(Fig.|7J. It is also possible that the source inflated bubbles of rising 
gas, carrying energy and pressure, as observed for example in M87 
(e.g Ghizzardi et al. 2004). However, radio sources are short-lived 
(10^ — 10* yr) and when they are starved of fuel they lose power. If 
this is the scenario we are observing, the radio source should be rel- 
atively old, and is in its fading phase. Hints in favour of this picture 
are the small size of the radio source, the distribution of the hard 
emission in the central 5 arcsec (Fig. 0, the small entropy core, 
and the relatively large external pressure at the location of the ra- 
dio galaxy. We cannot constrain the spatial distribution of the hard 
emission, which is strongly limited by the statistics. It is unlikely 
that the hard emission is non-thermal, since high energy electrons 
cannot live longer than their parent population. Thus it is possible 
that what we are observing is the signature of gas heating. 

The second possibility is that the source is young. The engine 
has just turned on and the source is expanding. This would explain 
i) the low surface brightness of the gas at the location of the radio 
source, and ii) the high surface brightness just outside it. This can 
be interpreted as an expanding source which is pushing the ambient 
gas away. If the source is young then it may not have had time to 
affect the entropy of the gas, which could explain the relatively 
small heating signatures in the entropy profile. 



7 CONCLUSIONS 

The Chandra observation of EXO 0423.4-0840 shows that a BL 
Lac object is absent and we observe an extended, cluster-scale 
source, which was previously discussed by Edge & Stewart (1991). 
The central X-ray emission is likely to be due to the ISM of the 
early-type galaxy at the centre of the cluster, and only at r > 40 
kpc does the ICM dominate over the ISM. The central X-ray emis- 
sion was little affected by heating in the way observed for other 
low-temperature clusters. If the source is relatively young, our re- 
sults represent indirect support for models where AGN activity in 
the central galaxy may be cause of the breakdown of self-similarity. 
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